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Abstract 

We offer a plausible resolution of the paradox (first formulated by Jevicki and Ro- 
driguez in Phys. Lett. B 146, 55 (1984)) that the two shifted harmonic oscillator 
potentials V{q) = + G/q^ + const may, in spite of their exact solvability in a 
non-empty interval of the couplings G, become supersymmetric partners if and only 
if G vanishes. We show that and how their G 7^ SUSY may be re-established via 
a regularization provided by pseudo-Hermitian quantum mechanics. 
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1 Motivation 



A suitable algebraic background of the theoretical construction of multiplets which 
would unify some experimentally observable bosons with fermions is provided by the 
graded Lie algebras, so called superalgebras. In such a setting, an exceptional role 
is played by the linear harmonic oscillator in one dimension, D = 1. Indeed, its 
Hamiltonian may be factorized and, subsequently, shifted to the left or to the right, 

ffiLHO) =p^^q2 = A- B- l = B- A + l, A = q + ip, B = q - ip 



^(LHO) _i=B-A, 



HiLHO) + i = A-B. 



According to Witten one may then introduce a "super-Hamiltonian" and two 
"supercharges" , 
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and notice that they generate a representation of Lie superalgebra sl(l/l), 

{Q, Q} = n, {Q, Q} = {Q, Q} = 0, [n, Q] = [H, Q] = 0. 

In the language of physics, one can speak about the bosonic and fermionic vacuum 
annihilated by both the supercharges. 



(g|0,0) 



exp(-gV2)/v/^ 




Q|0,0) = Q|0,0) = 0. 



The "bosons" themselves may be then introduced as created and/or annihilated by 
the first-order differential operators a'^ ~ 5 and/or a ~ A, respectively. The parallel 
creation, annihilation and occupation-number operators are also easily defined for 
"fermions" , 
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The supercharges become factorized as well, Q ~ ajF^ and Q, ~ a^jF. In terms of 
the harmonic-oscillator eigenstates \n) the Fock space will be spanned by the states 
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\nb-,nf) characterized by the presence of bosons and nj fermions where nj is equal 
to zero or one, 

|n,0) = 

In this way the harmonic oscillator may be understood as a next-to-trivial superym- 
metric field theory in a one- dimensional space-time which unifies the bosons with 
fermions. More details may be found, e.g., in the concise review paper 

In what follows, we shall analyze what happens if we replace the one-dimensional 
j^{LHO) 2;a^fjia[^ D— dimensional generalization with any real D > 2, 

= + + g^ a = iD-2)/2 + i, £ = 0,1,.... 

We shall be guided by the Witten's supersymmetric quantum mechanics where the 
use of the operators A = dg + W and B = —dg + W with arbitrary superpotentials W 
leads to the same supersymmetric pattern as above, forming the superHamiltonian 
from the two partner operators 

H^L) = B ■ A=p^ + - W, H^E) = A- B =p'^ + + W'. 

In the spirit of our recent letter we shall admit that these operators are pseudo- 
Hermitian |Q. 
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2 Bound states in the pseudo-Hermitian setting 

Beyond the elementary harmonic oscillator let us now contemplate a generalized 
superpotential 

W^"'\r) = r — '^^^ — — r = r(x) = x — ie. (1) 
r 

In the other words, we assume that we start from the choice of a real parameter 
7 and define the pair H{^l^r) of non-Hermitian operators. Such a recipe gives the 
partner Hamiltonians in the above D— dimensional harmonic oscillator form where 
Dl 1^ Dr in general, 

= ^^"^ - 27 - 2, i/S) = iJ^'^) - 27, a=|7|, P=\l + l\. (2) 
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In the light of ref. |^ the complexified line of coordinates r = x — i e circumvents 
the singularity in the origin so that the bound state wavefunctions are regular and 
expressible in terms of the Laguerre polynomials, 

Together with their energies 

E = E^^^ = AN + 2Q + 2, g = -Q.a, Q = ±1, N = 0,1,... 

these states are numbered by the integer N and by the so called quasi-parity Q. 



3 Supersymmetry, pseudo- Hermit ian way 

For reasons explained in ref. we must assume that 7 7^ 0,1,2,.... Up to that 
constraint, we may visualize the above construction as one of the regularizations 
recommended in the recent literature 0. Here we intend to summarize and discuss 
the subject in more detail. 

In the first step we notice that the quasi-parity Q coincides with the ordinary 
spatial parity P in the limit a 1/2. In such a limit the basis states are well known 
(cf. Appendix A). Once we move to a 7^ 1/2 we notice that the quasi-even states 
'ip(r) ~ r^/'^-"' still lie below their quasi-odd complements 'ijj(r) ~ 7-1/2+° fixed 
N. 

Whenever we choose a > 1, the limiting transition e — >■ moves the quasi- 
even solutions out of the Hilbert space completely. Otherwise, these states remain 
normalizable in a way depicted in Figure 1 where the following ordering is obtained 
for the A^— th bunch of the energy levels, 

e}-)") a{N)] < e'^^ h{N)] < E^^;^ c{N)] < ^f+f d{N)]. (3) 

This ordering is preserved along all their 7— dependent variation. Each of these four 
curves is just a once broken straight line but, in our picture, our eyes are guided 
by an infinitesimal shift of the levels in such a way that their shape may be easily 
followed (one should only recollect that the system remains undefined at all the 
integers 7 G IV). 
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In the Figure the physical, Hermitian hmiting transition e — has been per- 
formed. The general, e 7^ has been discussed elsewhere 0. We may only note here 
that in contrast to the latter and manifestly non-Hermitian, e 7^ scheme of ref. 
0, all our present states belong to the Hilbert space of the ordinary quantum me- 
chanics. Thus, our new scheme may be interpreted as a result of a pseudo-Hermitian 
regularization recipe studied, in more detail, elsewhere |^ (cf. also Appendix B for 
some more details). 

The inspection of Figure 1 reveals a certain generalized supersymmetry (SUSY) 
where the standard requirements of quadratic integrability tolerate the quasi-odd 
levels at all 7 but confine the existence of the levels a{n) to the very short interval 
of 7 G (—1, 1) and the existence of the levels b{n) to the interval of 7 G (—2, 0). As a 
consequence, one has to distinguish between the following five mutually significantly 
different regimes. 

1. "Far left" with 7 G {—00, —2) and with the complete degeneracy 

c{N)] = ) d{N)] 

where the Witten's index vanishes and where SUSY itself is broken because 
the ground state energy remains positive. All the spectrum is equidistant. 

2. "Near left" with 7 G (—2, —1) and with the mere partial degeneracy which 
survives from the preceding interval. There emerges the new series of ener- 
gies Ej^R)'^ [= b{N)] without any left partners; this possibility represents just a 
weaker form (and/or a more singular analogue) of the Jevicki- Rodriguez break- 
down of SUSY as mentioned above in Abstract. 

3. "Central domain" with 7 G ( — 1,0). This is the most interesting domain 
where the properties of the well known linear special case (which has 
7 = —1/2) appear generalized to the whole neighboring interval. Up to the 
exceptional (and newly emerging) ground state a(0) we may spot here the well 
known pattern of degeneracy, 

eI~-^ a{N)] < i^f^f ) b{N)] = E[lf cm < 
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<E'^!;> rf(iV)] = a(iV + l)<... 

so that SUSY becomes unbroken even for the spectrum which ceased to be 
equidistant at 7 7^ —1/2. 

4. "Near right" with 7 G (0, 1) and with the properties which "mirror" the far 
left (the series of energies Ejj^^^ [= b{N)] ceases to exist, etc). 

5. "Far right" with 7 G (1, 00), degeneracy 

and with the characteristic 7— independence of the almost completely degen- 
erate unbroken SUSY spectrum. 

We may summarize that the resulting SUSY pattern is fairly unusual. It may be char- 
acterized by several above-listed appeahng properties but one should re-emphasize, 
first of all, that near 7 = —1/2 a nice non-equidistant generalization of the textbook 
D = 1 SUSY oscillators is obtained. 
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Figure captions 

Figure 1. SUSY and the 7— dependence of the spectrum which is generated by the 
superpotential (p. 

Figure 2. Graphical solution of the selfconsistency condition E{q) = g in the 
schematic example of Appendix B. 3. with QN-k = 1 and fk = 3. 

(A) curve (p!TD at = 0.7 in the Hermitian regime. 

(B) curve (|ll|) at = 0.7 in the non- Hermitian regime (both energies are real), 

(C) curve ([TT|) at = 1.4 in the non-Hermitian regime (no real root, both 
energies are complex), 

(D) selfconsistency line E{g) = g. 
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Appendix A: The standard harmonic osciUator ba- 
sis on L2{M) 

Eigenstates of a Hamiltonian H{g) which commutes with the parity V may be num- 
bered by an integer n and by the superscript which characterizes the spatial parity 
of the state, 

Hig)\n^^\g))^Ei^\g)\n^^\g)) (4) 
These eigenstates form a complete basis in the Hilbcrt space L2(iR), 

oo 

E E \m^'^\g)){m^'^\g)\ = I . (5) 
The usual condition of their orthonormalization reads 

/oo 
{n^^\g)\x){x\m^'^\g))dx^SmnSa,r. 
-oo 

For the particular and exceptional harmonic oscillator H(0) = H^^^^^ these eigen- 
states are proportional to the well known Hermite polynomials, 

{x\n^+\Q)) = N2n'H2n{x) exp(-xV2) = 
(a;|n(-)(0)) = M2n+i'H2n+i{x) exp(-xV2) = (g) 



X e 




n = 0, 1,... . 



At each particular subscript n — m the pairs of the latter harmonic-oscillator basis 
states have an opposite parity, V\s) = +\s), V \t) = —\t). They may be transformed 
into the unitarily equivalent pairs of states 

( IT) ) ^ 72 ( i 1 ) ( It) ) 

with the real VT parities, VT\S) = +\S), VT \T) = —\T) where the complex 
conjugation T defined by the simple rule TiT = —i mimics the usual antilinear time 
reversal. 
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Appendix B: Main differences between the Hermi- 
tian and non-Hermitian Hamiltonians 



B. 1: P— symmetric models and the bases on L2(iR"^) 

Any eigenstate oi H = iJ^ = VHV satisfies its Schrodingcr equation even after 
a pre-multiplication by the parity V. Both the old and new eigenstates belong to 
the same real eigenvalue E which cannot be degenerate due to the Sturm-Liouville 
oscillation theorems. One of the superpositions 1-0) ± V\il)) must vanish while the 
other one acquires a definite parity. This is the essence of the mathematical proof that 
the V symmetry of wave functions cannot be spontaneously broken, V\n^'^\g)) — 



The knowledge of the spatial symmetry of the Hermitian Hamiltonian H[g) en- 
ables us to simplify many considerations and calculations by choosing and fixing the 
parity of the solutions in advance, 



This permits us to live, conveniently, on the semi-axis of x G (0, oo) = IR^ . In such a 
setting we rarely imagine that we are tacitly changing the Hilbert space from L2{M) 
to L2{]R^). We feel that this is a technicality which deserves a separate remark. 

On the new space (or, if you wish, in the old space equipped by the projector 
or singular metric IT) the inner product changes its meaning since we have to inte- 
grate over the mere semi-axis (symbolically, {'ip\'ip') — >■ {ip\Ii\ip')). This re-scales the 
orthogonality relations. 



Alternatively, we may omit the symbols H and switch to the two re-normalized bases 

(,x|(T„) = M2n'H2n{x) exp(-a;^/2), 

{x\Tn) = M2n+lH2n+\{x) exp(-xV2), 



±|nW(^)). 



{x\n^^\g)) = ±{{-x)\n^^\g)). 




ran- 
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which are both orthonormal on the half-hne. In parallel, condition splits in the 
two independent completeness relations 

oo oo 
m=0 n=0 

The overlaps of the states with different superscripts do not vanish and form a unitary 
matrix which changes the basis in L2{M~^). Its elements 

roc 

Un,m = 2 {n^+\^)\x){x\m^~\Q)) dx = 2{Sn\Ii\tm) = MTrn) (7) 







may be computed by the direct symbolic integration in MAPLE giving the exact 
values sampled in Table 1 from which one may extract some closed formulae, e.g.. 



(2n)!V%rf2 \ln + \T{n + \] 



22"+i(n!)2v^ TT r(n + l)' 
This means that the original vectors @ form an over-complete set and we may 
make a choice between the two alternative basis sets {|(t)} and {|t)}. They are both 
complete on the new Hilbert space L2{1R^). 

B. 2: PT— symmetric models 

The above-mentioned rigidity of the conservation of parity is lost during the transi- 
tion to the VT symmetric models H = = VTHVT where any quantity exp 
is an admissible eigenvalue of the operator VT since its component T is defined as 
anti-linear, Ti = —i. In more detail, every rule VT\tp) = exp{iip) lip) implies that 
we have 

VrVT\^) = exp{-iip)VT\tlj) = 

as required. The Schur's lemma ceases to be applicable. In the basis of Appendix A 
with the properties VT\S) = \S) and VT\L) = —\L), the general expansion formula 

oo 

H ^ ^ ^ \Sm}'^ m,n{Sn\ ~l~ \Lrri)Q'm,n{Ln\ ~l~ ^ \ Sm}Cm,n{Ln\ ~\~ ^ | -^m)^r?i,n ('S'n | ^ 

m,n=0 

contains four separate complex matrices of coefficients. Once it is subdued to the 
requirement H = VTHVT, we get the necessary and sufficient condition demanding 
that all the above matrix elements of H = must be real. 
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As long as the similar trick has led to the superselection rules for the spatially sym- 
metric Hamiltonians, we may conclude that the VT symmetric analogue of the direct- 
sum decompositions and superselection rule is just the much weaker constraint (§). 



B. 3: Schematic finite-dimensional matrix model with and 
without VT symmetry 







(9) 



Let us contemplate the partitioned matrix Schrodinger equation 

F-EI aA \ ( u 
At G-EI j ' \ w 

where F = F'^agger, G = and either a = 1 (Hermitian case) or a = — 1 {VT 
symmetric case). Schrodinger equations with the matrix representation generalize 
the models with VT symmetry [T^. Their spectrum may happen to be real and 
discrete, at least in the limit A —>■ 0, ot containing the complex conjugate pairs. Let 
us now descibe their nontrivial, non-perturbative solvable example. 

Preliminarily, both the Hermitian submatrices F and G of the Hamiltonian should 
be diagonalized via a pair of some suitable unitary transformations, F ^ f , G g. 
Their respective spectra {/„} and {gn} will be assumed real and discrete. 

Secondly, we shall ignore all the small elements of the coupling matrix A in our 
pre-diagonahzed effective Schrodinger equation (^, 

1 



f -EI-aA- 



■A^ 



y 



0. 



\ 9 -EI J 

Here, a = ±1 "remembers" its respective Hermitian and non- Hermitian origin and all 
the small elements in A are irrelevant causing just a perturbative, small deformation 
of the decoupled spectrum {/„} + {gn}- 

Thirdly, let us choose the latter coupling matrix in the form dominated by the 
transposed diagonal. 



/ 



A 









ai 





y Oat 







ao 







10 



Quite unexpectedly, this choice makes the problem exactly solvable since the secular 
equation det{Heff{g) — EI) = can be immediately factorized, 

0=f/o-^-a^Vf/i-^-«^^V------f/iv-^-«-''"''^ 



9N-QJ V 9N-1-QJ V 9o-Qj 

The explicit evaluation of zeros of the k—th factor is trivial, 

Eig)^fk-aak a^. (11) 

9N-k - Q 

The implementation of the selfconsistency g — E{g) gives the sequence of the mere 
quadratic algebraic equations 

1 

9N-k - E 



fk-E-auk al = 0. 



All their roots are available in closed form, 

Ek± = ^ (/fc + 9N-k ± \/(/fc-5'iv-fc)^ + 4Q;|afc|2^ . 

This confirms our a priori expectations since the Hermitian energies with a — 1 are 
always real while, at a — —1, we get the real spectrum if and only if 

\ak\ < \fk-9N-k\/2, A; = 0,l,...,iV. (12) 

Vice versa, we get a complex conjugate pair Ek± whenever we move to the strongly 
non-Hermitian regime and encounter a large and strong off-diagonality or coupling 
of modes in A. This is an independent linear-algebraic re-confirmation of the similar 
observations made during the explicit computations using differential equations. 
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Table 1. Overlaps (s„|n| t^) defined by eq. (0). Rows are numbered by 
n = 0, 1, . . . , 6, columns by m = 0, 1, . . . , 4. 



1/2^2 


-1/6^3 


1/20 




1/48^35 


1/2 


1/4 ^pi^f2 


-1/24 ^f^^/2 






-1/24 ^/2^fl 


1/8 VsVe 


1/16 v^Ve 


-^Vrov^ 




1/40 ^/2^/l 


-1/24 ^pi^fl 


1/16 v^VS 


1/32 VTOVS 


-^x/35v^ 
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